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ABSTRACT 

We present a new measurement of the mass-concentration relation and the stellar-to-halo mass 
ratio over a 5 x IO'^Mq to 2 x IO^^Mq mass range. To achieve this, we use the CFHT Stripe 
82 Survey (CS82) weak lensing data combined with a well defined catalog of clusters (the 
redMaPPer catalogue) and the LOWZ/CMASS galaxies of the Sloan Digital Sky Survey-Ill 
Baryon Oscillation Spectroscopic Survey Tenth Data Release (SDSS-III BOSS DRIO). The 
stacked lensing signals around these samples are modeled as a sum of contributions from 
the central galaxy, the dark matter halo, and the neighboring halos. We measure the mass- 
concentration relation: C 2 oo{M) - with A = 5.25 + 1.67, Z? = -0.13 + 0.12 for 

0.2 < z < 0.4 and A = 6.77 + 1.13, B = -0.15 + 0.06 for 0.4 < z < 0.6. We conclude that 
the amplitude A and slope B are both consistent with the simulation predictions by Klypin et 
al. (2014) within the errors. We also measure the stellar-to-halo mass ratio and find it to be 
flatter than previous measurement for high stellar masses because of the complex structures 
and merger history in massive dark matter halos. 

Key words: large-scale structure of Universe-gravitational lensing: weak-galaxies: clusters 
cosmology: theory dark matter 


1 INTRODUCTION 

Dark Matter (DM) is the dominant mass component in the universe. 
The hierarchical cold DM model with a cosmological constant 
(ACDM) is successful in explaining many observations on large- 
scale structures. The Navarro-Frenk-White (NFW) density profile 
(Navarro et al. 1997), which has been derived from DM only simu¬ 
lations, provides an accurate description of the dark matter density 
profile over most radii, can be described by just two parameters: 


* E-mail: huanyuan.shan@epfl.ch 


the concentration and the mass enclosed in a region with a mean 
density above some given threshold. 

Simulations have shown that these two parameters are cor¬ 
related, with the average concentration of a halo being a weakly 
decreasing function of mass (e.g. NFW; Bullock et al. 2001; Hen- 
nawi et al. 2007; Duffy et al. 2008; Prada et al. 2012). With the 
latest MultiDark simulations, Klypin et al. (2014) provide a more 
accurate analytical fit for dark matter halo density profiles and con¬ 
centrations. 

However, the mass-concentration (m-c hereafter) relation is 
different in the exact form as measured in simulations and esti- 
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mated from observations (e.g. Bullock et al. 2001; Comerford & 
Natarajan 2007; Hennawi et al. 2007; Duffy et al. 2008; Maccio et 
al. 2008; Mandelbaum et al. 2008; Oguri et al. 2009, 2012; Okabe 
et al. 2010; Sereno & Covone 2013; Foex et al. 2014). For example, 
considering a sample of 28 galaxy clusters whose radial profiles 
are well constrained from combining weak and strong lensing data, 
Oguri et al. (2012) found a steeper function of mass, with the slope 
-0.59 ±0.12. Merten et al. (2014) andUmetsu et al. (2014) stud¬ 
ied the m-c relation with the Cluster Lensing and Supernova Sur¬ 
vey with Hubble (CLASH) data, which is in excellent agreement 
with simulations when the CLASH X-ray selection function and 
projection effects are taken into account. Du & Fan (2014) studied 
the bias from the miscentering problem, selection effects and shape 
noise from intrinsic ellipticity of background galaxies on the mea¬ 
surement of m-c relation with individual clusters. They found that 
the existence of noise can compensate the miscentering problem 
and reduce the systematic bias although the scatters of mass and 
concentration get considerably larger. 

Mandelbaum et al. (2008) performed a statistical analysis to 
estimate the m-c relation over three orders of magnitude in mass 
using the stacked weak lensing signals of galaxies, groups and clus¬ 
ters in Sloan Digital Sky Survey (SDSS). Their analysis used the 
information outside the halo central region in order to avoid re¬ 
gions where baryons dominate, and to minimize errors from mis- 
estimation of clusters’ centers. They found an m-c relation with a 
slope consistent with simulations, but an amplitude that is about 2cr 
below theoretical expectations. With the SDSS-III/Constant Mass 
(CMASS) galaxies in the Canada-France-Hawaii Telescope Lens¬ 
ing Survey (CFHTLenS), Miyatake et al. (2013) also showed a 
20% smaller concentration than the theoretical predictions. With 
the CFHTLenS catalog and photometrically selected clusters, Cov¬ 
one et al. (2014) found an over-concentrated but flat relation. 

Galaxy properties, such as luminosity or stellar mass, are 
tightly coupled to the DM halo mass. The link between the stellar 
and DM halos is critical for understanding the formation and evo¬ 
lution of galaxies. The observed luminosity of the central galaxy 
is always used to fix the stellar mass component. There are several 
methods to measure the stellar-to-halo mass relation (SHMR) in¬ 
cluding stellar kinematics (Conroy et al. 2007; More et al. 2011), 
abundance matching (Moster et al. 2010; Behroozi et al. 2010), the 
Tully-Fisher relation (Pizagno et al. 2007; Miller et al. 2014) and 
weak lensing (Mandelbaum et al. 2006; Leauthaud et al. 2012; Ve- 
lander et al. 2014; Hudson et al. 2015; Han et al. 2015). In particu¬ 
lar, Leauthaud et al. (2012) found that the slope of the SHMR rises 
sharply at M, >~ 4.5 x 10'°Mq. 

In this paper, we use a sample of redMaPPer clusters (Rykoff 
et al. 2014), the LOWZ and CMASS galaxies of SDSS-III BOSS 
DRIO and the shear catalog from the CFHT Stripe 82 dataset to 
measure the tangential shear, and quantify the m-c relation and the 
SHMR. 

2 OBSERVATIONAL DATA 
2.1 The Source catalog 

The source galaxies used in our measurement are taken from the 
CFHT Stripe 82 Survey (CS82: Erben et al. 2015), which is an i- 
band imaging survey covering a large fraction of SDSS Stripe 82 
region with a median seeing 0.59”. The survey contains a total of 
173 tiles with 165 from CS82 and 8 f rom CFHT-LS Wide. Each 
CS82 tile was obtained from 4 consecutive dithered observations 


each with an exposure time of 410 seconds. The completeness mag¬ 
nitude is Iab ~ 24.0. 

The shapes of faint galaxies are measured with the lensfit 
method (Miller et al 2007, 2013). In this work we use the same 
pipeline as the CFHTLenS collaboration. Therefore the shear cali¬ 
bration factors described in Heymans et al. (2012) can be directly 
applied to our data. 

We use all galaxies with weight w > 0, FITCLASS=0 
and MASKED in which w represents an inverse variance weight 
assigned to each source galaxy by lensfit, FITSCLASS is a 
star/galaxy classification provided by lensfit. The parameter MASK 
describes the mask information at an object’s position. Objects with 
MASK<1 can safely be used for most weak lensing analyses (Er¬ 
ben et al. 2013). After masking out bright stars and other image 
artifacts, the effective sky coverage is ~130deg^. 

As CS82 contains only /-band, we derive the photometric 
redshifts (photo-z) using EAZY (Brammer et al. 2008) for the 
source galaxies from the overlapping multi-color co-add data from 
SDSS (Annis et al. 2014). The total number of source galaxies 
that pass the aforementioned cuts and have an available photo-z 
is 4,449,618 (~ 9.Sgals jarcmin^). 


2.2 Lens selection 

We use the redMaPPer clusters and LOWZ/CMASS galaxies in two 
redshift bins, 0.2 < z < 0.4 and 0.4 < z < 0.6, as lenses: 

(1) The RedMaPPer Cluster Catalog: 

The redMaPPer cluster catalog by Rykoff et al. (2014) is based 
on the optimized red-sequence richness estimator A. The richness 
A is the number of red sequence galaxies brighter than 0.2L, at 
the redshift of the cluster within a scaled aperture. The redMaPPer 
cluster catalog in CS82, which is obtained from the SDSS Stripe 82 
coadd using the CS82 object detections, covers the redshift range 
0.1 < z < 0.6. The objects with richness T > 15 correspond to the 
halo mass ~ 8 x 10'^ Mq (Erben et al. 2015). In order to reduce the 
effects of miscentering problem, we choose the clusters with well- 
defined centers Pcen > 0.9, which is the center probability. We then 
divide the clusters into two redshift bins with redshift 0.2 < z < 0.4 
and 0.4 < z < 0.6. For both low and high redshift clusters, we 
divide them into two richness bins (20 < T < 30,30 < T < 200) 
with (57,31) and (96,75) systems, respectively. 

(2) LOWZ & CMASS Galaxy Catalog: 

We also use the LOWZ and CMASS galaxy catalog from 
SDSS-III BOSS DRIO as lens galaxie^ The two samples are 
both the main SDSS-III BOSS BAO tracers (Dawson et al. 2013). 
LOWZ samples consist of red galaxies at z < 0.4 from the 
SDSS DR8 (Aihara et al. 2011) imaging data; CMASS samples 
are selected with an approximately constant stellar mass thresh¬ 
old (Eisenstein et al. 2011). Stellar masses of the two samples 
are from the Portsmouth SED-fitting (Maraston et al. 2013) based 
on the BOSS spectroscopic redshift. We choose the LOWZ and 
CMASS galaxies within redshift 0.2 < z < 0.4 and 0.4 < z < 0.6. 
We also divide them into 3 and 4 stellar mass bins, which will 
be used on the SHMR measurement, with (2350,2492,2325) and 
(5966,6168,7481,7611) systems for the low-z and high-z range, 
respectively (see Eigure 1). 


* https: //WWW. sdss3.org/drlS/spectro/galaxy.php 
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3 LENSING SIGNAL 


We measure the m-c relation with stacked leasing signal. Stack¬ 
ing many halos reduces the fluctuations due to noise caused by 
uncorrelated structures along the line-of-sight, shape measurement 
noise, substructures, and the shape variations of individual halos. 
Thus, the measurement can determine the average mass profile. 
Furthermore, it allows for the leasing measurement of lower mass 
halos, where individual detection is impossible due to their smaller 
shears relative to clusters. Individual massive cluster observations 
and those based on stacked analysis of many halos are thus com¬ 
plementary, drastically increasing the available baseline in mass. 

Except the dark matter halos, the models for the lensing sig¬ 
nals include the following components: 

(1) Gavazzi et al. (2007) found that a central baryonic compo¬ 
nent is required to fit both the strong and weak lensing profiles of 
early-type galaxies in the SLAGS survey. This component can be 
simply modeled as a point mass (e.g., Leauthaud et al. 2012). 

(2) In order to measure the lensing signals, locating the centers 
of dark matter halos is very important because miscentering can 
add significant systematic uncertainties. In general, halo centers are 
typically assumed to be traced by a central galaxy. However, it has 
been argued that these assumptions may not be entirely true and 
that central galaxies could be “sloshing” around in the dark matter 
halo potential well (e.g. van den Bosch et al. 2005; Gao & White 
2006; Johnston et al. 2007; Hikage et al. 2013). As in Johnston et 
al. (2007), Leauthaud et al. (2010), George et al. (2012) and Li et 
al. (2014), centroid errors will lead to a smoothing of the lensing 
signal on small scales and to an underestimate of halo masses. The 
distribution of miscentering can be described by a 2D Gaussian 
distribution: 


P(Roff) = — exp[-^(J?off/cro(f)A (1) 

o-off \ 2 I 

where Raff is the distance between the halo center and central 
galaxy position, (Toff is the effective scale length. This model is 
from mock catalogs by Johnston et al. (2007), which is consistent 
with the measurement from 2D lensing shear maps of 25 massive 
clusters by Oguri et al. (2010). 

(3) At the large scales, the neighboring halos contribute to the 
lensing profile (Johnston et al. 2007; Cacciato et al. 2009, Li et al. 
2009, Oguri & Hamana 2011). The tangential shear profile due to 
the neighboring halos is 


y2h{0\M,z)= -hild)- - 

J 2n (1 -Hz)3EcritD2(z) 


( 2 ) 


where Ji is the second order Bessel function, Pm(z) is the mass 
density at z, D/^{z) is the angular diameter distance, F^ffc) is the 
linear matter power spectrum, Xcrit is the critical mass density, and 
bk{M) is the halo bias that we take from Tinker et al. (2010). 

Combining all the contributions, the lensing signal AS can be 
written as: 


AS(J?) - + PcenASNFwJJ^) + (3) 

(1 - Pcen)AENFw(^l^off) + 


where sequentially, the terms come from the contribution of central 
baryonic components, the dark matter halos, the miscentering ha¬ 
los, and the neighboring halos. The model is described by 5 param¬ 
eters: the halo mass Myir, the concentration Cvk, the central baryonic 
mass Mp, the effective scale length (Xofr and the probability pcen that 
the centers we use here are accurately centered on the dark matter 
halos. 


To obtain AS we stack lens-source pairs in 8 logarithmic radial 
(R) bins from 0.1 to 10 Mpc h“*. For a sample of selected lenses, 
A1,(R) is estimated using 


AZiR) = 


2uls ^hyt ^crit 
2/1 Wls 


(4) 


where yj' is the tangential shear, W/., = and w„ is a weight 

factor introduced to account for intrinsic scatter in ellipticity and 
shape measurement error. 

With the multiplicative errors m taken into account statistically 
(Miller et al. 2013), we apply it to the average shear measurement in 
both the low-z and high-z samples using the correction as Velander 
et al. (2014) and Hudson et al. (2015): 


1 + K(z,) = 


YjhWi.d +m) 


The lensing signal can be calibrated as follows: 


AZ‘^\R) = 


AZ(R) 

l+R(z,)' 


(5) 

(6) 


The effect of this correction term is to increase the average lensing 
signals by ~ 5.5% and 6.2% at 0.2 < z < 0.4 and 0.4 < z < 0.6, 
respectively. 

To minimize the contamination due to photo-z uncertainties of 
the source galaxies, we remove the lens-source pairs with Zs - Zi < 
0 . 1 . 


4 RESULTS 

We show the observed galaxy-galaxy lensing signals for the low 
and high redshift samples in Figure 1. The black circles with er¬ 
ror bars represent the measured excess surface mass density. We fit 
them to the model described in section 3. The black solid curves 
shows the best-fit results with all the components. The different 
model components are shown in dashed curves: the central bary¬ 
onic matter (red), the dark matter halo profile (green), the miscen¬ 
tering halo component (blue) and neighboring halos (magenta). At 
small scales (0.1 - IMpc), the dark matter halo is the dominant 
component. At large scales (> 3Mpc for galaxies and > 5Mpc 
for clusters), the neighboring halos contribute more. The best-fit 
results are given in Table 1. We convert the measured Myi, and 
Cvir to M 200 and C200 with the formula by Johnston et al. (2007). 
The stellar masses M, of LOWZ/CMASS galaxies are from the 
Portsmouth SED-fitting (Maraston et al. 2013) based on the SDSS- 
III BOSS DRIO data. M, of the redMaPPer clusters are estimated 
by matching the positions of central galaxies with the galaxy cat¬ 
alog of SDSS DRIO. The fitted probability pcen are all close to 1, 
which indicates that: (1) the centers of redMaPPer clusters are well 
defined; (2) LOWZ/CMASS galaxies are almost all central galax¬ 
ies, ~ 10% of them may be the most luminous satellite galaxies 
as the clustering analysis by Parejko et al. (2013) and White et al. 
(2011). 

4.1 m-c relation 

The m-c relation from the simulation predictions is given as: 

C 2 oo(z,Af)=A|^J (l-rz)'", (7) 

with Mo = 2.0 X lO'^Mg/t^*. As we cannot get a redshift depen¬ 
dence relation with only two redshift bins, the redshift dependence 
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0.2<z<0.4 


0.4<z<0.6 




0.1 1 0.1 1 0.1 1 0.1 1 0.1 1 

R [Mpc/h] R [Mpc/h] 


Figure 1. Best-fit model for the low (left) and high (right) redshift samples (black solid lines). Black circles show the measured excess surface mass density 
AE of both LOWZ/CMASS galaxies and clusters with different bins (clusters: Bottom-middle and bottom-right panels; LOWZ/CMASS: others). The model 
components are the central baryonic matter (red), the dark matter halo profile (green), the miscentering halo component (blue) and neighboring halos (magenta). 


Table 1. Density Profile Models in Figure 1. 



Z 

Mvir 

^^vir 

M 200 

C 200 

M, 

<^ofF 

Pccn 



I0‘^Mo//i 


10'*Mo/h 


lO'^Mo 

Mpc 


(10*'“, 10"-^)Mo 

0.2 < 2 < 0.4 

0.12 ±0.04 

4.66 ± 1.29 

0.10 ±0.03 

3.80 ± 1.09 

0.13 ±0.10 

0.45 ± 0.26 

0.76 ±0.18 

(J0l1.2,i0ll.35)^^ 

0.2 < ^ < 0.4 

0.15 ±0.04 

4.19 ±0.85 

0.13 ±0.03 

3.41 ±0.71 

0.19±0.17 

0.30 ±0.17 

0.82 ±0.13 

(10ii.35_ ioi2.0)Me 

0.2 < ^ < 0.4 

0.37 ± 0.07 

3.60 ±0.66 

0.31 ±0.06 

2.91 ± 0.55 

0.30 ± 0.27 

0.19±0.15 

0.89 ±0.10 

20 < 4 < 30 

0.2 < z < 0.4 

1.01 ±0.09 

3.15 ± 1.02 

0.83 ± 0.073 

2.54 ± 0.85 

0.32 ± 0.25 

0.55 ± 0.28 

0.95 ±0.12 

30 < i < 200 

0.2 < z < 0.4 

1.37 ±0.36 

3.66 ± 1.20 

1.14 ±0.30 

2.96 ± 1.01 

0.37 ± 0.21 

0.63 ± 0.25 

0.96 ±0.11 

(10*““, 10“‘)Mo 

0.4 < z < 0.6 

0.074 ± 0.03 

5.40 ± 1.51 

0.067 ± 0.03 

4.62 ± 1.31 

0.095 ± 0.056 

0.32 ±0.17 

0.86 ±0.15 

(10***, 1O“-2)M0 

0.4 < z < 0.6 

0.083 ± 0.03 

5.00 ± 0.97 

0.074 ± 0.03 

4.26 ± 0.84 

0.14±0.11 

0.23 ±0.12 

0.87 ±0.13 

(10‘*^10“-3^)Mo 

0.4 < z < 0.6 

0.12 ±0.05 

4.69 ±0.91 

0.11 ±0.04 

3.99 ± 0.79 

0.19 ±0.15 

0.40 ±0.17 

0.86 ±0.12 

(I0‘'-3^10‘2-“)Mo 

0.4 < z < 0.6 

0.24 ± 0.06 

3.88 ±0.82 

0.22 ± 0.04 

3.28 ± 0.70 

0.30 ±0.18 

0.24 ±0.14 

0.88 ±0.21 

20 < 4 < 30 

0.4 < 2 < 0.6 

0.68 ±0.11 

3.62± 1.11 

0.60 ± 0.084 

3.06 ± 0.96 

0.33 ± 0.21 

0.57 ± 0.33 

0.91 ±0.21 

30 < 4 < 200 

0.4 < z < 0.6 

1.39 ±0.30 

3.35 ±0.52 

1.20 ±0.26 

2.83 ± 0.45 

0.47 ± 0.28 

0.68 ± 0.21 

0.90 ± 0.33 


Notes: M, is the stellar mass from the Portsmouth SED-fitting (Maraston et al. 2013) based on the SDSS-III BOSS DRIO data. 


is calculated with the mass and concentration with Multidark sim¬ 
ulations of Klypin et al. (2014). Because of the upturn of the con¬ 
centrations at large masses regime of high redshift, the redshift de¬ 
pendence C-0.67 is fitted with the DM halos of Multidark in the 

redshift region 0.0 < z < 0.6, which is suitable to the data in our 
analysis. The best-fit m-c relations are shown in Table 2. Because of 
the bigger errors on low-z sample, the m-c relations of both samples 
are consistent. 

Giocoli et al. (2012) found that the halo triaxiality and the 
substructures within the host halo virial radius can bias the m-c 
relation. They propose a method for correcting the m-c relation for 
the projection etfects and also for adiabatic contraction: 

C2d(M) = C3d(M) X 1.630M-'’ '’‘^ (8) 

In Figure 2, we show the results on the 3D corrected m-c relation 
after rescaling the results to z = 0. The 3D correction is used on all 
the observational dataset and simulation predictions. 


Table 2. Best-fit Mass-Concentration relations. 



z 

A 

B 

0.2 < 

z < 0.4 

5.25 ± 1.67 

-0.13 ±0.12 

0.4 < 

z < 0.6 

6.77 ± 1.13 

-0.15 ±0.06 


With the 3D corrections by Giocoli et al. (2012), we can di¬ 
rectly compare our measurements with the three-dimensional (3D) 
N-body simulation results. Comparing with the results by Duffy et 
al. (2008) and Klypin et al. (2014), the amplitudes A and slopes B 
of both redshift samples are consistent with the simulation predic¬ 
tions within the errors. 

We then compare our results with previous measurements. The 
3D corrections by Giocoli et al. (2012) are used on these measure¬ 
ments. 

By stacking weak lensing signals of galaxies, groups and clus¬ 
ters from SDSS dataset, Mandelbaum et al. (2008) fit the data to the 
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Figure 2. 3D corrected m-c relation with clusters and galaxies in CS82 after 
rescaling all the results to j = 0. Red and blue circles denote the stacked 
lensing signals from low and high redshift samples, respectively. The red 
curve is the best-fit m-c relation of the low redshift sample. The blue curve 
and cyan-shaded area is best-fit m-c relation of the high redshift sample and 
its Icr uncertainty. Black curves are the simulation predictions by Duffy 
et al. (2008) (dashed) and Klypin et al. (2014) (solid). The black symbols 
denote the measurement on mass and concentration by Mandelbaum et al. 
(2008) with SDSS; Covone et al. (2014) with CFHTLenS; Umetsu et al. 
(2014) and Merten et al. (2014) with CLASH cluster sample. 


Figure 3. Halo-to-stellar mass ratio as a function of stellar mass. Red and 
blue dots denote the measurements from low and high redshift samples, re¬ 
spectively. The red solid curve and cyan-shaded area is the best-fit SHMR 
relation of our low redshift sample combined with HI5 and its Icr uncer¬ 
tainty. The blue solid curve is the best-fit SHMR relation of our high redshift 
sample combined with H15. The red and blue dashed curves are the fitted 
relation by the COSMOS measurements We compare our measurements 
with other galaxy-galaxy lensing measurements, including SDSS (Mandel¬ 
baum et al. 2006), CFHTLenS (H15), COSMOS (Leauthaud et al. 2012), 
COSMOS X-ray clusters (Leauthaud et al. 2010), and GAMA (Han et al. 
2015). 


model assuming a spherical NFW profile excluding small scales to 
reduce the effects of baryons and miscentering halos. As shown 
in Figure 2, their measurements are comparable to ours, but more 
scattered. 

With photo-z selected clusters (Wen et al. 2009) in the 
CFHTLenS data, Covone et al. (2014) measured the mass and con¬ 
centration by stacking shear profiles. The measured concentrations 
in four richness bins are higher than the simulation predictions and 
generate a flat and over-concentrated mass-concentration relation. 

A possible reason might be the choice of the cluster catalog. 
Wen’s cluster finder purely based on the photo-z of galaxies. The 
richness of some clusters may be overestimated due to the inaccu¬ 
racy of photo-z. A detail comparison will need to be carried on in 
future. 

The recent CLASH measurements for individual clusters 
(Merten et al. 2014) and stacking (Umetsu et al. 2014) focus on 
more massive strong lensing clusters ~ IO'^Mq/j^*, which give 
large concentrations. Accounting for the CLASH selection func¬ 
tion based on X-ray morphology and projections effects inherent in 
lensing observations (Meneghetti et al. 2014), a better agreement 
can be found between observations and theoretical predictions. 


4.2 SHMR measurement 

With the measurement of stellar mass, we can also estimate the 
SHMR. Following Leauthaud et al. (2012), fsHMR(Mh) is mathe¬ 


matically defined via its inverse function: 

loglo(/5^«MK(^A)) = logloW) +/?logl0 + (9) 


where Mi is a characteristic DM halo mass, M,fi is a characteristic 
stellar mass, j8 is the low-mass SHMR slope, S and y is the high- 
mass SHMR slope. In Figure 3, we show the results of this paper in 
addition to other galaxy-galaxy lensing studies, including the mea¬ 
surement in SDSS (Mandelbaum et al. 2006), CFHTLenS (Hud¬ 
son et al. 2015 (H15)), the Cosmic Evolution Survey (COSMOS: 
Leauthaud et al. 2010, 2012) and the Galaxy And Mass Assembly 
(GAMA: Han et al. 2015). Our measurement of M 200 IM, is ideally 
combined with the results of H15, as they provide very good esti¬ 
mates of the SHMR at the same redshift bins as ours and covering 
smaller masses. Our best fit results for this combined dataset are 
given in Table 3. 

We compare our results with the best-fit in Leauthaud et al. 
(2012). At low stellar mass {M, < 4 x 10'°Mo), the constraints 
are only from H15. The behaviors of SHMR are similar. The major 
halo mergers are minor galaxy mergers in this regime. The growth 
of stellar mass through the effects of mergering will be limited com¬ 
pared to the growth of DM halo mass. At M, ~ 4 - 6 x IO'^Mq, 
M 200 IM, of the results reaches a minimum, which reflects the most 
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Table 3. Best-fit results of SHMR. 


z 

log 10 (Ml) 

logio(M.,o) 

ys 

5 

r 

0.2 < z < 0.4 

12.51 ±0.047 

10.97 ± 0.032 

0.47 ±0.019 

0.43 ±0.11 

1.45 ±0.27 

0.4 < z < 0.6 

12.69 ± 0.056 

11.13 ±0.035 

0.48 ± 0.022 

0.44 ±0.15 

1.90 ±0.31 


efficient star formation rate at such a mass. The minimum values 
are consistent with the COSMOS measurements within the errors. 
At low stellar mass regime, the DM halos are dominated by single 
central galaxies, the supernovae feedback could be the main rea¬ 
son to prevent the star formation. At high stellar mass, we find a 
relatively flatter slope. The scatter should be due to the appearance 
of hot gas in the massive DM halos, making more complex merger 
history. 


5 CONCLUSIONS 

Combing two different redshift samples of redMaPPer clusters and 
LOWZ/CMASS galaxies with the shear catalog of CS82, we fit 
the lensing signals around these clusters and massive galaxies in¬ 
cluding the effects of baryons, miscentering halos, and neighboring 
halos, and study the m-c relation and also the SHMR. 

Our main results are as follows: 

1. We find that the concentration decreases weakly with mass, 
c = A(M/Mo)® with the amplitude A and slope B in good agree¬ 
ment with the previous measurement by Mandelbaum et al. (2008). 
Including the 3D correction suggested by Giocoli et al. (2012), 
the two parameters are also consistent with simulation predictions 
within the errors. 

Klypin et al. (2014) also shows that the NFW profile can sys¬ 
tematically overpredict the halo concentration by 10 - 20% for the 
massive halo density profile. The Einasto profile can provide a more 
accurate measurement. However, the current precision of data can¬ 
not distinguish between the two profiles. The lensing signals in the 
inner region of dark matter halos have large error bars. Mandel¬ 
baum et al. (2008) found a very similar m-c relation for the NFW 
and Einasto profiles with SDSS data. Using the recent CLASH, 
Umetsu et al. (2014) found similar mass and concentration of the 
two profiles. In order to compare the two profiles, we also fit the 
lensing signals of the redMaPPer clusters with the Einasto profile, 
and find similar m-c relation as with the NFW profile. We conclude 
that the current data can not distinguish the two profiles. 

2. Combining previous galaxy-galaxy lensing measurements, 
we measure a relatively flatter slope of the SHMR at high stellar 
mass. With COSMOS data, Leauthaud et al. (2012) found a steeper 
slope of the SHMR at M, > 4.5 x 10 '®Mq. 

The flatter slope of the SHMR at high stellar mass implies the 
massive DM halos contain more complex structures and merger 
history. The scatter at the high stellar mass regime could be due to 
the appearance of hot gas in the massive DM halos, making a more 
complex merger history. 

In our analysis, we use the stellar mass of central galaxy 
for the redMaPPer clusters. There could be more than one central 
galaxy in a cluster, however, making our measurement a lower limit 
on the total stellar mass in a halo. Such an effect could raise the 
SHMR at high stellar mass regime and increase the scatter. 

With future wide surveys (DES/HSC/LSST/EUCLID), we ex¬ 
pect to measure more accurate the m-c relation and the SHMR, and 
test for possible measured bias. 
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